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ABSTRACT 

Recent observations of Type la supernovae (SNe la) suggest that some of the progenitor white dwarfs (WDs) 
had masses up to 2.4-2.8 M Q , highly exceeding the Chandrasekhar mass limit. We present a new single 
degenerate (SD) model for SN la progenitors, in which the WD mass possibly reaches 2.3-2.7 M Q . Three 
binary evolution processes are incorporated; optically thick winds from mass-accreting WDs, mass-stripping 
from the binary companion star by the WD winds, and WDs being supported by differential rotation. The WD 
mass can increase by accretion up to 2.3 (2.7) M Q from the initial value of 1.1 (1.2) M Q , being consistent 
with high luminosity SNe la such as SN 2003fg, SN 2006gz, SN 2007if, and SN 2009dc. There are three 
characteristic mass ranges of exploding WDs. In an extreme massive case, differentially rotating WDs explode 
as an SN la soon after the WD mass exceeds 2.4 M because of a secular instability at r/|W| ~ 0.14. For 
a mid mass range of Mwd = 1.5-2.4 M Q , it takes some time (spinning-down time) until carbon is ignited to 
induce an SN la explosion after the WD mass has reached maximum, because it needs a loss or redistribution of 
angular momentum. For a lower mass case of rigidly rotating WDs, M WD = 1.38-1.5 M , the spinning-down 
time depends on the timescale of angular momentum loss from the WD. The difference in the spinning-down 
time may produce the "prompt" and "tardy" components. We also suggest the very bright super-Chandrasekhar 
mass SNe la are born in a low metallicity environment. 

Subject headings: binaries: close — stars: winds, outflows — supernovae: individual (SN 2003fg, SN 2007if, 
SN 2009dc) — white dwarfs 



1. INTRODUCTION 

Type la supernovae (SNe la) play important roles in astro- 
physics, e.g., standard candles in cosmological distance and 
main production cites of iron-group elements. However, the 
nature of SN la progenito rs is not fully understood y et (e.g., 
iNiemever & HiUebrandtll2004t iNomoto et al. 1120001) . It has 
been widely accepted that SNe la are explosions of carbon- 
oxygen (C+O) white dwarfs (WDs). The observed features of 
SNe la are better explained by the Chandrasekhar mass model 
rather than the sub-Chandrasekhar mass model. It is not clear, 
however, how the WD mass gets close enough to the Chan- 
drasekhar mas£] for carbon ignition; i.e., whether the WD 
accretes H/He-rich matter from its binary companion [single 
degenerate (SD) scenario], or two C+O WDs merge [double 
degenerate (DD) scenario]. 

Recent observations of several very bright SNe la suggest 
that their progenitor WDs might have super-Chandrasekhar 
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1 Hereafter we us e the Chandrasekhar mass Mch = 1.46(Y C /O.5) 2 M0 (e.g. 
[Chandrasekhar 19393 for a non-rotating WD with electron mole number F c - 
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Super-Chandrasekhar mass of the p rogen itors of SNe la 
was first reported by iHowell et all ((2006) on SN 2003fg 
(~ 2.1 M C+O WD with 1.29 ±0.07 M Q 56 Ni). More 
candidates of super-Chandrasekhar mass progenitors were 
added to t he list, i.e., SN 2006gz, SN 2007if, and 
SN 2009dc. iHillebrandt et al. I d2007l) suggested a possibil- 
ity that SN 2003fg is an aspherical explosion of a Chan- 
drasekhar mass WD. However, based on their recent spec- 
tropolarimetric observation of SN 2009dc which sho ws lit- 
tle continuum polarization of < 0.3 Tan aka et al. I (120101) con- 
cluded that a very large asphericity was not likely to occur 
at the explosion of SN 2009dc. Moreover, the very bright 
SN 2007if require a 2.4 ± 0.2 M prog enitor C+O WD with 
1.6 ±0.1 M Q 56 Ni dScalzo et al. 1120101) and SN 2009dc also 
deman ds > 2.0 M Q progeni tor C+O WD with 1.4-1.7 M Q 
56 Ni dSilvermanetl^l20Til) or ~ 2 8 M Q C+O WD with 
~ 1.8 Mq 56 Ni dTaubenberger et al. I I20TTI) . Such a large 
WD mass, highly exceeding Mch, could challenge both the 
DD and SD scenarios. Some authors suggest a DD merger 
scenario for these super-Chandrasekhar mass SNe la (e.g., 
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IScalzo et al. 1 120 lot ISilverman et al. 112-01 ll) . mainly because 
the existing SD scenario provides only M W d ^1-8 M Q (e.g., 
iChen & Li ll2009tlLiu et al. l20Toh . 

As for theoretical ex plosion models of super- 
Chandrasekhar mass WDs, iSteinmetz et al. I (119921) per- 
formed simulations of pure detonation models and found that 
almost all the matter is burned into iron group elements. Their 
results could not explain the observed interme diate mass ele- 
ments seen in the spectra of SNe la. Recently. iPfannes et al. I 
(2010) revised the simulations of pure detonation models 
based on Yoon & Langer's (2005) super-Chandrasekhar mass 
models. They found that not all the WD matter is burned into 
iron group elements and some intermediate mass elements 
are synthesized in the outer low density region of the WD. 
The synthesized mass of nickel is large enough (~ 1.5 Mq) 
to explain the super luminous SNe la such as SN 2009dc, 
which supports the super-Chandrasekhar mass models. 

In the present paper, we show that such high WD masses 
are possibly explained by the SD scenario. There are three 
key binary evolution process es, i.e., optical l y thic k winds 
from mass-accreting WDs (Ha chisu et al. I 1 1996b . mass- 
strippi ng from the bin ary companion star by the WD winds 
(Hach isu et al. Il2"008al) . and differential rotation which makes 
possible a WD more m assive than the Chandrasekhar limit 
(lYoon & Langer] l200l . The main difference of our SD 
model from the existing SD scenarios is that we take into 
account the effect of mass-stripping by the winds from mass- 
accreting WDs. We describe our basic assumptions and meth- 
ods of binary evolutions in Section 2. Our numerical results 
are presented in Section 3. A possible binary evolution to our 
super-Chandrasekhar mass SNe la is suggested in Section 4. 
Discussion and concluding remarks follow in Sections 5 and 
6, respectively. 

2. MODELLING OF MASS-STRIPPING EFFECT AND BINARY 
EVOLUTION 

First we briefly describe the method of our calculations 
for binary evolutions which is essentially the same as in 
lHachisu et al. I (l2008al) except that the evolution of the com- 
panion star is computed by a Henyey-type code. In a binary, 
the more massive (primary) component evolves to a red gi- 
ant star (with a helium core) or an asymptotic giant branch 
(AGB) star (with a C+O core) and fills its Roche lobe. Dy- 
namically unstable mass transfer begins from the primary to 
the secondary and forms a common envelope. After this first 
common envelope evolution, the separation shrinks and the 
primary becomes a helium star or a C+O WD. The helium 
star evolves to a C+O WD after a large part of helium is ex- 
hausted by core- and shell-helium-burning. Then the binary 
consists of a close pair of a C+O WD and a main-sequence 
(MS) star; from such a pair we start our calculations. 

When the secondary expands to fill its Roche lobe, mass 
transfer begins from the secondary to the WD. If the sec- 
ondary mass is more massive than the WD (precisely, when 
the mass ratio q = Mix* /M\$ = Mms.o /Mwd,o > 0.8), the mass 
transfer takes place on a thermal timescale. When the mass 
transfer rate exceeds the critical rate, 

M cl - = 6.68 x 10 -7 ( ~T7~~ ~ 0-445 J M© yr" 1 , (1) 

optically thic k winds blow from the WD (see lHachisu et al. I 
[1996111 999allbl: iNomoto et al. 1120071 for more details). 

The winds from the W D collide with the secon dary and 
strips off its surface layer (Hachisu & Kato 2003a b). This 



mass-stripping effect attenuates the mass transfer from the 
secondary to the WD, thus preventing the formation of a com- 
mon envelope even for a rather massive secondary. Thus the 
mass-stripping e ffect widens the donor mass range of SN la 
progenitors (e.g., [Hachisu et al. 2008aJ3). 

We hav e incorporated t he mas s-stripping effect in the same 
way as in lHachisu et al. I d2008al) ; i.e., the mass stripping rate 
M st ,ip is related to the wind mass-loss rate M w i n d as 

Mstrip=ClM„ind, (2) 

with the factor c\ defined as 

= ffeff • g(q) ( v 2 a \ n 
C1 "^l3-^ms \2GMj' W 

where M = M W d+M M s, M W d is the WD mass, M M s is the 
main-sequence companion mass, a is the separation of the bi- 
nary; (f>]_i and </>ms denote the Roche potential (normalized 
by GM/a) at the MS surface and the L3 point near the MS 
companion, respectively, i.e., </>L3 - </>ms ~ 0.3; 7j e ff is the ef- 
ficiency of conversion from kinetic energy to thermal energy 
by the shock, g(q) ~ 0.025 is the geometrical factor of the 
MS surface hit by the wind including the inclination (oblique 
shock) effe ct of the wind velocity against the companion's 
surface (see Hachisu et al. 1999a, for more details on g(q)), 
and q = Mz/M\ = M ms /Mwd is the mass ra tio. Here we mod- 
ify equation (21) of lHachisu et al. I (Il999al) to include the ef- 
fect of Roche lobe overflow from the L3 point. We assume 
?7eff = 1 in the present calculation. If the wind velocity is 
as fast as 4,000 km s~' we have ci ~ 10 as estimated by 
lHachisu & Katol d2003ab . lHachisu & Katol (2003a b) found 
best fit models at a = 1.5- 10 for RX J0513.9-6953 and at 
c\ =7-8 for V Sge, respectively. Here we adopt c\ = 1, 3, 
and 10 to check the dependence of mass-stripping effect on 
the parameter c\. 

The winds depend on the WD mass and metallicity (chem- 
ical composition). We have revised metallic ity effects of WD 
winds (Figure 1 of Kobavash i et al. |[l998l) . The results are 
shown in Figure Q] We did not find wind solutions but in- 
stead static envelope solutions for the "no winds" region. On 
the other hand, the "strong winds" region means that the 
win d velocity exceeds the e scape velocity (see, e.g., Figure 
4 of iKato & Hachisu"! 1 19941 for wind solutions). The veloc- 
ity of the "weak winds" does not exceeds the escape veloc- 
ity. Winds are accelerated mainly by the peak of OPAL opac- 
ity due to iron lines while they are sometimes accelerated by 
a small opacity peak due to helium ionization for the rela- 
tively low metallicity region. Here we regard that stripping 
effect works only for strong winds. We have strong wind so- 
lutions (i.e., stripping effect) for Mwd ^ 0.6 M© in the solar 
metallicity environment while for Mwd 1-0 Mq in a low 
metallicity region of Z = 0.001 (Population II), and further 
for Mwd ^1-2 Mq in the much lower metallicity region of 
[Fe/H]<-2. 

Mass accretion onto the WD spins it up because angular 
momentum is added to the WD (e.g . . lLanger et al. 1 120001; 
PaczvnskTI 1 1 99 It iPopham & Naravanl fl 99 1 ; Piersant i et al. I 
2003albt lUenishi et al. II2003I) . If the WD rotates rigidly, the 
mass of the WD can only slightly exceed the Chandrasekhar 
mass limit of no rotation, Mqh- More massi ve WDs seem t o 
be possible if they rotate differentially (e.g., Hachisu 1986). 
Based on simplified numerical calculations, lYoon & Langer 
(2004) showed that the WD can increase its mass beyond Mch 
when the accretion rate onto a WD is high enough (Mwd ^ 
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Fig. 1 . — Wind region in the metallicity-WD mass diagram. Data are essentially the same as those in Figure 1 of Kobayashi et al. ( 1998) but a few points are 
revised. Each symbol shows as follows: strong winds (blue circles), weak winds (green triangles), both of which are driven by the OPAL peak at logT (K) = 5.2, 
weak winds (green squares) which are driven by the opacity peak due to helium ionization, and no winds (red crosses). "Strong winds" mean that the wind 
velocity exceeds the escape velocity while "weak winds" do not exceed. The blue dotted line shows a border of strong winds. The red dash-dotted line indicates 
a border of no winds. In the low metallicity region, strong winds blow only for very massive WDs (e.g., Mwd > 1.0 Mq for [Fe/H]< -1 and Mwd > 1.2 Mq 
for [Fe/H]< -2). 

10~ 7 M Q yr -1 ) and the WD is differentially rotating. They 
showed that the gradient of rotational rate is kept around the 
critical value for the dynamical shear instability, because the 
timescale for the Eddington-Sweet meridional circulation is 
too long (~ 10 9 yr) to redistribute angular momentum in the 
WD core, and that the differential rotation is strong enough 
for the WD mass to exceed Mch significantly. 

In the present study, we simply assume that the WD is sup- 
ported by differential rotation and its mass can increase with- 
out carbon being fused at the center as long as the mass ac- 
cretion rate is higher than 3 x 1O~ 7 M yr _1 (= Mb). We have 
adopted the limit by considering the fact that when Mwd < Mb 
hydrogen shell-burning occurs intermittently and recurrent 
nova out bursts eject a large part of the hydrogen-rich enve- 
lope (e.g. Hach isu & Kato 11200 ll) . As a result, the net growth 
rate of the C+O core mass is reduced significantly so that the 
timescale of the angular momentum deposition on the WD 
would become comparable to the timescale for meridional cir- 
culation to re-distribute the angular momentum. 

3. NUMERICAL RESULTS 

3.1. Evolutions of WD+MS Systems 

We have investigated evolutions of binaries which ini- 
tially consist of a WD (primary) and a zero-age main- 
sequence (ZAMS) star, by calculating the evolut ions of the 
secon dary stars using a Henyey-type code (e.g., iSaio et ai~l 
1988). We have adopted a standard chemical composi- 
tion of (X,Z) = (0 70,0 2) and used OPAL opacity tables 
(Igl esias & Rogers | [l996). The mixing length is assumed to 
be 1 .5 times the pressure scale height. 

As the star evolves, the radius increases and eventually fills 



its critical Roche lobe. Then, a mass transfer from the sec- 
ondary to the primary WD begins. We have determined the 
rate of mass loss from the secondary, M%, at each stage by 
requiring for the stellar radius to be equal to the mean radius 
of the Roche lobe. If \M^\ is smaller than the critical rate 
M cr given in Equation ([T), all the matter from the secondary 
will be accreted onto the WD; i.e., M\ = \M2\- When the mass- 
accretion rate exceeds M a -, an optically thick wind blows from 
the WD. It causes mass-stripping from the secondary star as 
discussed in the previous section. Then, the systemic mass 
loss at a rate of M w j n d + M str j p occurs, which affects the evo- 
lution of the Roche lobe radius and hence M^. During the 
evolution, Mi is giv en as min(M cr , {Mil) (see lHac hisu et al. 
fl99HrT999lf2008al for details). 

In the previous work ([Hachisu et al. 2008a), we stopped 
the binary evolution at the epoch when the WD mass reached 
1.38 Mq, because we reg arded that a WD explodes as an 
SN la at M W D = 1 -38 M Q (lNomotolll982al) . In the present 
work, however, we further follow the binary evolution un- 
til the mass accretion rate of the primary WD decreases to 
Mi = M b = 3 x 10" 7 Mq yr" 1 , because the WD is expected to 
be supported by a strong differential rotation until then. 

Although Equation O is based on the assumption of spher- 
ical symmetry dNomoto et al. 1120071) . we consider it approxi- 
mately applicable even for Mwd > 1 .38 M Q for the following 
reasons. Equation (Q~|) is closely related to the local Eddington 
luminosity (and nuclear energy conversion rate per unit mass 
for hydrogen mass fraction, X = 0.70), which is mainly deter- 
mined by the local opacity (peak) at a part with temperature 
of ~ 2 x 10 5 K in the envelope. Since this layer is far from the 
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FIG. 2. — Binary evolutions of two WD+MS pairs: time stalls when the mass transfer begins from the MS to the WD. (a) A WD of Mwd.o = 1.1 Mq and a MS 
star of Mms.o = 4.7 Mq with the initial orbital period of P or b o = 0.543 day (near ZAMS). The red thick solid line: the mass-loss rate from the MS companion, 
— M%. The red thin solid line: one tenth of the mass-loss rate, — 0. IM2. The other solid lines denote the companion mass (blue), M2 = Mms, the orbital period 
(green), P or b, and the WD (primary) mass (black), Mi = Mwd, each symbol of which is attached to each line. The red dash-dotted line: stripping rate from the MS 
companion, M stl j p . The long red dashed line: the wind mass-loss rate from the WD, M wmL \. The red dotted line: the mass-increasing rate of the WD, M] = Mwd, 
which is merged to the red solid line in the later time. The mass-loss rate, — M2, increases up to 6 X 10~ 6 Mq yr - ' and quickly drops to 1 X 10~ 6 Mq yr -1 . 
Then, the optically thick wind stops and the mass-loss rate decreases to —Mi = 3 X 10~ 7 Mq yr -1 at t = 1.7 X 10 6 yr. At this time, the WD mass has increased 
to Mi = M W d = 2.3 Mq. (b) Same as (a) but for a pair of 1.2 Mq WD and 5.5 Mq MS with the initial orbital period of P orb ,o = 0.568 day (near ZAMS). The 
mass-loss rate increases up to — M2 = 1.5 X 10~ 5 Mq yr -1 and then quickly drops to 1 X 10~ 6 Mq yr -1 . Then, the optically thick wind stops and the mass-loss 
rate decreases to -M 2 = 3 X 10~ 7 Mq yr~' at t = 1 .8 X 10 6 yr. At this time, the WD mass has increased to M\ = Mwd = 2.7 Mq . 



C+O core surface (at the hydrogen shell burning region), the 
local Eddington luminosity should be insensitive to a rapid 
rotation of the C+O core. 

We have followed binary evolutions for various sets of the 
initial parameters (Mwd.o, M2,o,f i-b,o) and obtained maximum 
WD masses. In the first example, we set c\ = 10 and assume 
Mwd.o =1.1 Mq, M2.0 = 4.7 Mq, and P O rb,0 = 

0.54 day (the 

companion fills its Roche lobe after it slightly evolved off 
from ZAMS). The result is shown in Figure |2^. The mass- 
loss rate of the companion quickly increases up to — M2 = 
8 x 10" 6 M© yr" 1 and then drops to -M 2 = 1 x 10" 6 M Q yr" 1 
at t = 10 6 yr (— M2 is also plotted in the one tenth scale in the 
same figure by a red thin solid line). In the early phase of 
mass transfer, we have -0. IM2 ~ M^md because —M<i w M str j p 
and ci = 10. As mass transfer rate decreases, the optically 
thick wind stops at t ~ 1 x 10 6 yr, and the accretion rate M\ 
decreases to 3 x 10" 7 Mq yr" 1 at t = 1.7 x 10 6 yr. The WD 
mass has increased to 2.3 Mq at this epoch. (Even if we adopt 
Mb = 1, 2, 3, and 4 x 10~ 7 Mq yr" 1 , the maximum WD mass 
does not change so much, i.e., 2.45, 2.36, 2.29, and 2.22 Mq, 



respectively. The dependence of the WD mass on the value of 
Mb is very weak so that the adopted value itself is not essential 
for the results.) 

In the second example, we assume a more massive bi- 
nary with ci = 10, M WD .o =1.2 Mq, M 2 .o = 5.5 Mq, and 
^orb,o = 0.57 day (the companion fills its Roche lobe near 
ZAMS as in the first case). The result is shown in Figure 
|2j5. The mass-loss rate of the companion increases up to 
—Mi = 1.4 x 10~ 5 M Q yr" 1 and then quickly decreases. The 
optically thick wind stops at almost the same time as in 
Figure |2^, t = 1 x 10 6 yr, and afterward the mass transfer 
rate gradually decreased to M = M, = 3 x 10" 7 M Q yr" 1 at 
t = 1 .8 x 10 6 yr. The WD mass has increased to 2.7 M Q at this 
epoch. 

The third and forth examples are for the binaries in which 
the companion fills its Roche lobe in a later phase, i.e., near 
the end of central hydrogen burning (Figure HJ. The third 
example is for c\ = 10, M W d.o = 0.9 M Q , M 2 .o = 3.5 M Q , 
and P or b.o = 1 -38 day. The result is shown in Figure [3^. 
The mass-loss rate of the companion increases up to — M2 = 
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FIG. 3. — Same as Figur^2] but for the companion models near the end of central hydrogen-burning, (a) A WD of Mwd.o = 0.9 Mq and a MS of Mms.o = 3.5 Mq 
with the initial orbital period of P or t,.o = 1-38 day. The mass-loss rate from the companion decreases to — M2 = 3 X 1CT 7 Mq yr~' at ; = 1.67 X 10 6 yr. At this 
time, the WD mass has increased to 1.8 Mq. (b) A WD of Mwd.o = 10 Mq and a MS of Mms.o = 4.0 Mq with the initial orbital period of P or (, = 1.44 day. The 
mass-loss rate decreases to -Mi = 3 X 10~ 7 Mq yr -1 at t = 1.59 X 10 6 yr. At this time, the WD mass has increased to 2.0 Mq. 



3.8 x 1O~ 6 M0 yr -1 and then quickly decreases. The optically 
thick wind stops at t = 9 x 10 s yr, and afterward the mass 
transfer rate is gradually decreasing to Mb = 3 x 10~ 7 Mq yr -1 
at t = 1.7 X 10 6 yr. The WD mass has increased to 1.8 M Q at 
this epoch. 

In the fourth example, we have adopted c\ = 10, Mwd,o = 
1.0 M , M 2 , = 4.0 M , and P orb , = 1.44 day. The result 
is shown in Figure ^p. The mass-loss rate of the compan- 
ion increases up to — Mj = 6 x 10~ 6 M Q yr" 1 and then quickly 
decreases. The optically thick wind stops at t = 8 x 10 5 yr, 
and afterward the mass transfer rate is gradually decreasing to 
M b = 3 x 10" 7 M Q yr -1 at t = 1 .6 x 10 6 yr. The WD mass has 
increased to 2.0 M Q at this epoch. 

In all the models, the mass-loss rate from the secondary, 
—Mi, increases to > 3 x 10~ 6 M Q yr" 1 . If all the hydrogen- 
rich matter were converted into helium and then into carbon- 
oxygen and therefore accreted at the same rate as — M2, carbon 
burning would ignite off-center and th e C+O core would b een 
converted into O+Ne+Mg core (e.g.. iKawai et al. Ill988h so 
that the WD would not explode as an SN la. However, the net 
accretion rate to the WD is given as |M 2 1 - |M str i p | - |M w ; nc j 
(=Mi = Mwd = the dotted lines in Figures [2] and 0, not ex- 
ceeding 3 x 10~ 6 Mq yr" 1 . So we expect no off-center burning 
of carbon during the accretion phase and that the C+O core 
grows until carbon ignites at the center. 



Our numerical results are summarized in Figured] This fig- 
ure shows the maximum WD mass MwD.max obtained in the 
binary evolution against the initial companion mass M^o for 
various sets of parameters, Mwd.o, ^orb.o, and c\. Figure |4^ 
shows examples for the parameter of c\ = 10. We set a step of 
0.5 M Q for the initial companion mass (except for a few mod- 
els). The rightmost ends of these solid lines correspond to the 
model beyond which the mass transfer rate of —Mi quickly 
increases to 10~ 4 -10~ 3 Mq yr -1 and thus a common envelope 
is formed. We regard that two stars are merging into one. 
Solid lines are for the case that the companion fills its Roche 
lobe near ZAMS, i.e., P or b ~ 0.5 day. The WD mass reaches 
MwD.max ~ 2.7, 2.3, 1.9, 1.7, and 1.5 Mq, for the initial WD 
mass of M WD ,o = 1.2, 1.1, 1.0, 0.9, and 0.8 M Q , respectively. 
The dashed lines show the case that the companion fills its 
Roche lobe near the end of central hydrogen burning, i.e., 
^orbo ~ 1-5 day. In this case, Mwd max ~ 2.3, 2.2, 2.0, 1.8, 
and' 1.6 Mq, for M WD , = 1.2, 1.1, 1.0, 0.9, and 0.8 Mq, re- 
spectively. For the latter case (dashed lines) MwD.max smaller 
than that for the former case (solid lines) even for the same c\ , 
Mwd,o, andM2.o. This is because the more evolved secondary 
has a larger mass transfer rate, -M 2 , and thus the duration of 
a high mass-accretion phase, in which period the WD grows 
quickly, is shorter than that of the near ZAMS case. 

Figure [4J} shows examples for the parameters of c\ = 3 
(solid lines) and c\ = 1 (dotted lines), in which the companion 
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FIG. 4. — The maximum WD mass against the companion mass for the various initial WD masses and orbital periods, (a) ci = 10. The solid lines with small 
filled circle: the companion is slightly evolved off from ZAMS but still very close to ZAMS (P rb.o ~ 0.5 day) when the companion starts mass-transfer. The 
dashed lines with small open circle: the companion is close to the end of central hydrogen-burning (P or b.o ~ 1.5 day) when the companion starts mass-transfer. 
The horizontal thin solid line denotes Mwd = 1-38 Mq. (b) Same as those in (a) but for C[ = 3 (solid lines) and C[ = 1 (dashed lines), where the companion fills 
its Roche lobe near ZAMS, i.e., P or b.o ~ 0-5 day. See text for more details. 



fills its Roche lobe near ZAMS, i.e., P 01 -b,o ~ 0.5 day. The ob- 
tained maximum W D mass is smaller than that for c\ = 10. In 
our previous work (Hachisu et al. 2008a), based on a sim- 
plified model of mass-transfer, we showed that more mas- 
sive companions can avoid merging for larger c\ -parameters. 
In the present work, we adopted a more realistic compan- 
ion model which results in more accurate mass-transfer rates. 
Comparing Figure with Figure |4j?, we confirm that more 
massive companion avoids merging for a larger c\ -parameter. 

The maximum WD mass depends strongly on c\ and 
Mwd,o- To summarize, the maximum WD masses are 
MwD^max = 1-9, 1.9, and 2.3 Mq for cj = 1, 3, and 10, if 
the binary evolution starts from Mwd.o = 1.1 M Q . If the bi- 
nary evolution starts from a more massive initial WD mass of 
Mwd,o = 1 .2 M©, a larger maximum WD mass is resulted as: 
M W D,max = 2.0, 2.1, and 2.7 M for c\ = 1, 3, and 10, respec- 
tively. 

3.2. Initial-final States of Binary Evolutions 

We plot, in Figure [5] the initial parameter region that pro- 
duces super-Chandrasekhar mass (Mwd > 2.0 Mq) SNe la 
in the logP-M2 (orbital period - secondary mass) diagram 
(hatched by magenta lines). Here we assume the metallicity of 
Z = 0.02 and the initial white dwarf mass of Mwd.o = 1.1 Mq . 
A binary system initially inside the hatched region (above) 
increases its WD mass up to a super-Chandrasekhar mass of 



Mwd > 2.0 Mq. The companion mass decreases with time, 
so that its position moves downward. When it reaches the 
hatched region (below), the mass-transfer rate drops below 
Mb = 3 x 10~ 7 M Q yr" 1 and the Eddington-Sweet meridional 
circulation would redistribute angular momentum in the WD 
to trigger central carbon burning. Evolutionary paths of two 
typical cases of M2.0 = 4.0 Mq companion star near ZAMS 
(Po = 0.53 day) and M2.0 = 4.5 Mq companion near the end 
of central hydrogen burning (Po = 1 -48 day) are followed by a 
thick magenta line with an arrow. Our calculated final states 
fall into a very narrow strip because all their evolutions stop 
at the same mass transfer rates of M, = 3 x 10~ 7 Mq yr" 1 . 

For comparison, we also added our previous results taken 
from lHachisu et al. I (l2008albL Il999al) . where we assumed 
that a non-rotati ng C+O WD explodes at M W d =1-38 M Q 
dNomotolll982al) . The initial binary system inside the region 
encircled by a thin solid line (labeled "initial") increases its 
WD mass up to the critical mass (Mwd = 1.38M Q ) and ex- 
plodes in the regions encircled by a thick solid line (labeled 
"final") for the white dwarf and main-sequence star (WD + 
MS) system (left) and the white dwarf and red giant (WD + 
RG) system (right). We have found no super-Chandrasekhar 
mass SN la (M W d > 2.0 Mq) region for the WD+RG sys- 
tems mainly because mass of the companion is too small to 
supply much matter to the WD. Currently known positions 
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FIG. 5. — The regions that produce super-Chandrasekhar mass (Mwd > 2.0 Mq) SNe la are plotted in the logP— M2 (orbital period - secondary mass) plane 
(hatched by magenta lines). Here we assume the metallicity of Z = 0.02 and the initial WD mass of Mwd.o = 1.1 Mq. The initial system inside the hatched 
region (above) increases its WD mass up to a super-Chandrasekhar mass of Mwd > 2.0 Mq and then reaches the hatched region (below) where the mass-transfer 
rate drops below M/, = 3 X 10~ 7 Mq yr . Evolutionary paths of two typical cases (4.0 Mq companion star near ZAMS and 4.5 Mq companion near th e end o f 
central hydrogen burning) are traced by a thick magenta line with an open arrow. We also added the previous results taken from Hachisu et al. 1 2008a b, 1999a), 
where we assumed that a non-rotating C+O WD explodes at Mwd = 1-38 Mq. The initial system inside the region encircled by a thin solid line (labeled "initial") 
increases its WD mass up to the critical mass (Mwd = 1 -38M0 ) and explodes in the regions encircled by a thick solid line (labeled "final") for the white dwarf and 
main-sequence star (WD + MS) system (left) and the white dwarf and red giant ( WD + RG) system (r ight). Currently known posi tions of Galactic recurrent nova e 
and supersoft X-ray sourc es are indicated by a star mark (★) for U Sco (e.g., [H achisu et al. 2000a), a triangle for T CrB I Belczyriski & Mikolaiewska HT998l) , 
a filled circle for RS Oph iBrandi et al. 2009), a square for V Sge (Hachisu & Kato 2003b). Arrows indicate V394 CrA and CI Aql, because their companion 
masses are not known. 



of some Galactic recurrent novae and supersof t X-ray sources 
are ind icated by a star mark (*) for U Sco ("e.g.JHachisu et alTI 
2000a), a triangle for T CrB (Belczyriski & Mikolaiewska 
19981) . a filled cir cle for RS Oph (IBrandi et al. 1 120091). a 
square for V Sge dHachisu&Katoll2003bl) . and by arrows 
for the two recurrent novae, V394 CrA and CI Aql, because 
mass of a companion is not yet accurately obtained. Two sub- 
classes of the recurrent novae, the U Sco type and the RS Oph 
type, correspond to the WD + MS channel and the WD + RG 
c hannel of SNe la, respec tively. 

lHachisu & Katol (200 1 ) summarized the mass transfer rates 
of Galactic recurrent novae in their Table 2, in which only U 
Sco has a large mass transfer rate enough to support differen- 
tial rotation, i.e., the mass transfer rate is ~ 2.5 x 1O~ 7 M0 yr" 1 
and the net growth rate of C+O core is ~ 1.0 x 10~ 7 M Q yr" 1 . 
There fore, U Sco may have a super-Chandrasekhar mass 
WD. iThoroughgood et al. I (1200 11) observationally estimated 
the WD mass of U Sco to be M WD = 1 .55 ± 0.24 M in the 
1999 outburst. The systematic error is, however, too large to 
draw a definite conclusion for super-Chandrasekhar mass, so 
we strongly encourage more accurate observations to deter- 
mine the WD mass in U Sco. 

3.3. Stability of Mass-accreting, Differentially Rotating WDs 

In the present paper, we simply assume that the mass- 
accreting, differentially rotating WDs are stable and can ac- 
crete mass from the companion as long as its mass accretion 
rate is larger than M\, = 3 x 10~ 7 Mp, yr" 1 . Th is assumption is 
based on the results in lYoon & Langerl ((2004). However, they 
did not incorporate the possibility of baroclinic instability in 



the WD mainly because equi-density and equi-pressure lines 
are essentially parallel in the WD du e to high degene racy. If 
the baroclinic instability could occur (Fujimoto 1993), angu- 
lar momentum is effectively transported in a shorter timescale 
than the accretion tim escale, so that a rigid rotation is realized 
dSaio & No moto 2004). If this is the case, suggested super- 
Chandrasekhar mass WDs would not be realized in the SD 



Recently, iPirol (120081) concluded, based on Fujimoto's 
(1993) baroclinic instability criterion, that baroclinic insta- 
bilities bring the WD very close to solid-body rotation dur- 
ing accretion phase, which makes it difficult to grow a WD 
mass substantially past the Chandrasekhar mass. However, 
we should notice that Fujimoto's criterion on baroclinic in- 
stability is just a necessary condition but not a sufficient con- 
dition. In this sense, Piro's (2008) result is not always cor- 
rect. We must keep in mind that he showed only a possibil- 
ity for baroclinic instability on mass-accreting WDs but not 
an occurrence of baroclinic instability. On the other hand, 
the Richardson criterion on the Kelvin-Helmholtz instability 
is also a necessary condition but has ever been tested for a 
variety of basic (plane-parallel) flows as a suffic ient condi- 
tion for the instability (see, e.g., Fuiimoto 1993). Yoon & 
Langer's (2005) stability criterion was based on the Richard- 
son criterion of the Kelvin-Helmholtz instability. Therefore, 
we expect that the true mixing instability for redistribution of 
angular momentum occurs somewhere between the Richard- 
son criterion and Fujimoto's (1993) criterion. 

If we assume Yoon & Langer's (2005) stability criterion 
(the Richardson criterion) and extrapolate their formula on 
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FIG. 6. — Schematic illustration for evolutions of the WD mass against time. Time starts when the secondary fills its Roche lobe and it is switched from linear 
to logarithmic at t Rd 3 X 10 6 yr (a vertical magenta line). The evolutional lines in the left (part) to the magenta line are taken from our numerical results while 
those in the right (part) are just schematic illustrations. From top to bottom: (a) Mwd.o = 1-2 Mq, Mi,q = 5.5 Mq, Pq = 0.57 day, ci = 10. Same model as the 
model in Figure|2p. In this model, T /\ W\ exceeds 0.14 before the mass accretion rate becomes lower than Mb = 3 X 10~ 7 Mq yr~' . We expect an SN la explosion 
at the position of magenta star mark (Mwd ~ 2.4 Mq) soon after a secular instability sets in. (b) Mwd.o = 1-1 Mq, Mi q = 4.7 Mq, Pq = 0.54 day, C[ = 10. 
Same model as the model in Figure|5Ji, where 7//|W| does not exceed 0.14 before M] < Mb = 3 X 10~ 7 Mq yr . Loss or redistribution of angular momentum 
leads to an increase in the central density until carbon is ignited to induce SN la explosion (Mwd ~ 2.3 Mq). The delay time depends on the timescales of 
angular momentum los s or redistribution. It may take a time of order 10 8 -10 5 yr if angular momentum redistribution occurs due to Eddington-Sweet circulation 
I Yoon & Langer 2004). (c) M WD , = 1.0 Mq, M 2 ,q = 4.5 Mq, Pq = 1.46 day, c, = 10. The WD mass reaches Mwd = 1-85 Mq, which is differential-rotation- 
supported, (d) Mwd.o = 0-9 Mq, M20 = 3.5 Mq, Pq = 1.38 day, C| = 10. The WD mass reaches Mwd = 1-76 Mq, which is differential-rotation-supported. 
(e) Mwd.o = 0-8 Mq, Mi q = 4.9 Mq, Pq = 0.52 day, c\ = 10. The WD mass reaches Mwd = 1-47 Mq, which can be supported by rigid rotation (Hachisu 
1986; Uenishi et al. 2003). Therefore, it may take a longer time than the case of differentially-rotating-supported WDs. (f) Mwd.o = 0.7 Mq, M20 = 3.0 Mq, 
Pq = 0.47 day, C] = 5. The WD mass reaches Mwd = 138 Mq, which can also be supported by rigid rotation. 



(their Equation [29]) to higher masses, where T and 
W are the rotational and gravitational energies of the WD, 
respectively, angular momentum redistribution could occur 
when T/\W\ exceeds 0.14 and a secular instabilit y sets in. 
This instability occurs at M WD ss 2.4 M Q (see also lHachisu I 
1986). At this time, the WD would be too massive to 
keep hydrostatic equilibrium and contracts on a much shorter 
timescale due to angular momentum transport. As a result 
of such a rapid compression, the WD eventually ignites ex- 
plosive carbon burning in the central region to give rise to an 
SN la explosion. If super-Chandrasekhar mass SNe la occur 
in this way, the WD mass can be as massive as up to 2.4 Mq. 
To confirm this prediction, accurate estimates of the ejecta 
masses, in which asphericity is taken into account, would be 
important. 

3.4. Evolutions in Three Mass Ranges of Exploding WDs 

If we apply the above discussion to the super- 
Chandrasekhar mass WDs, we may further discuss the 
fate of the WDs. There are three characteristic mass ranges 
of exploding WDs depending on the secular instability and 
rotation law: (1) Mwd > 2.4 M Q (a secular instability at 
T/\W\= 0.14), (2) M WD = 1. 5-2.4 M Q (differential rotation), 
and (3) M W d =1.38-1.5 M Q (rigid rotation). 

3.4.1. SNe la triggered by Secular Instability: 
SN 2007 if /SN 2009dc-like 

For very massive WDs, a secular instability sets in at 
T/\W\ = 0.14 and the WD explodes as an SN la when it 
reaches 2.4 M Q . For example, the model of Mwd.o = 1-2 M Q , 



M 2 fl = 5.5 Mq, P = 0.57 day, and c { = 10 (Figure|2p) reaches 
T/\W\ = 0.14 at f = 1.28 x 10 6 yr. This evolution model is also 
shown in Figure [6]by the uppermost line a ("a" is attached to 
the line). Therefore, we expect an SN la explosion to occur 
at the magenta star mark because the timescale of the secular 
instability is much shorter than the mass-accreting timescale. 
This kind of binary evolution corresponds to the very lumi- 
nous super-Chandrasekhar mass SNe la such as SN 2007if 
and SN 2009dc. It is interesting that these brightest super- 
Chandrasekhar mass SNe la are the youngest among our mod- 
els (younger than a few hundred Myr). 

3.4.2. SNe la triggered by Loss or Redistribution of Angular 
Momentum: SN 2003fg/SN 2006gz-like 

For M WD = 1.5-2.4 M Q , T/\W\ does not exceed 0.14. For 
example, the model for the initial parameters of Mwd,o 
I , I M Q , M 2 ,o = 4.7 Mq, P = 0.54 day, and cj = 10 (Figure 
|2k) reaches Mwd = 2.3 Mq before the accretion rate drops to 
Mi < M b = 3 x 10" 7 M Q yr" 1 at the red filled circle in Fig- 
ure [6] (line b). We also show two other examples. One is 
for the initial parameters of M W d,o =1-0 M Q , M 2 ,o = 4.5 M Q , 
Po = 1 .46 day, and c\ = 10. In this case, the WD mass reaches 
Mwd = 1 -85 M Q , which is also differential-rotation-supported 
(line c in Figure [6). The other is for Mwd,o = 0.9 M Q , 
M 2 , = 3.5 M , Pq = 1.38 day, and c y = 10. The WD mass 
reaches Mwd = 1 -76 Mq, which is again differential -rotation- 
supported (line d in Figure [6j. 

However, the WD does not explode soon after the WD mass 
has reached the maximum mass (1.5-2.3 Mq), because the 
central density is low so that central carbon is not ignited yet. 
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After some time, loss or redistribution of angular momentum 
would lead to an increase in the central density of the WD. 
As a result, carbon is ignited to induce an SN la. This wait- 
ing time depends entirely on the timescale of angular momen- 
tum loss or redistribution. It is an order of 10 s — 10 9 yr if an- 
gular momentum redistribution occurs d ue to the Eddington- 
Sweet circulation (lYoon & Langerll2004b . If direct loss of an- 
gular momentum occurs in a much shorter timescale (due to 
magneto-dipole radiation or so), it could be much shorter than 
~ 10 8 yr. In any case, it is shorter than 10 s — 10 9 yr. Therefore 
we suggest that these relatively young SNe la correspond to 
the "prompt" component. We regard that this kind of binary 
evolution corresponds to super-Chandrasekhar mass SNe la 
such as SN 2003fg/SN 2006gz, not so extremely luminous as 
SN 2007if/SN 2009dc. 

3.4.3. SNe la triggered by Loss of Angular Momentum: normal 

SNe la 

For Mwd < 1.5 Mq, WDs could be rigidly rotating. For 
example, line e in Figure|6]shows the case of Mwdo = 0.8 M Q , 
M 2 , = 4.9 M , Po = 0.52 day, and c\ = 10. The WD mass 
reaches Mwd = 1-47 which can be support ed by rigid 
rotation (e.g.. lHachisTnil986l:IUenishi et al.1l2003b . Also, line 
f in Figure [6] depicts the case of initial parameters of Mwd o = 
0.7 M Q , M 2 ,o = 3.0 M Q , P = 0.47 day, and c x = 5. The WD 
mass reaches Mwd = 1 -38 M Q at the red filled circle. The WD 
can also be supported by rigid rotation. 

The WD will eventually explode as an SN la after the WD 
spins down and the central density increases high enough to 
ignite carbon. This spinning-down timescale depends highly 
on the subtraction mechanism of angular momentum from the 
WD such as magnet o-dipole radiation in a strongly magne- 
tized WD. Recently, lllkov & Soker I (12011b reexamined the 
timescale of spin-down due to r-mode gravitational radiation 
and showed that r-mode i nstability is not significant in spin- 
ning down WDs (see also Lindblom 1999: lYoon & L anger 
2005). Instead, they suggested that magneto-dipole radiation 
leads to spinning-down in a typical timescale of < 10 9 yr in 
WDs as massive as Mwd ^ 1 -6 M Q but in a typical timescale 
of > 10 9 yr in M WD = 1.4-1.5 M Q . If it is the case, these 
rigidly -rotating WDs of Mwd = 1.38-1.5 M Q may correspond 
to the "tardy" component, suggesting that the companion has 
evolve d off to a white dwarf when the WD explodes as an 
SN la (Uusthamll20TTtlDi Stefano et al. Il20ll . 

There is an example of WDs that had once been extremely 
spun-up but now are very slowly spinning down. Recently a 
massive, fast spinning WD with a spin period of P sp j n = 13.2 s 
and a mass of Mwd = 1-28 ±0.05 Mq was found in a bi- 
nary system HD 49798/RX J0648.0-4418 (P orb = 1.55 day), 
which is consisting of a hot subdwarf and a WD (e.g., 
Mereg hetti et al. I [201 lb . This spin period is only a factor 
of two longer than the critical spin periocfl The spin-down 
rate is too small t o be d etected, i.e., f sp i n < 9 x 10"° s s" 1 . 
Mereghetti et al71 (l201 lb estimated a very weak magnetic field 
of B < 1 kG at the WD surface. In such a low magnetic 
field, the spin-down rate due to magneto-dipole radiation 
should be so small that we expect a very long decay time 

2 The total angular momentum and rotational energy are J = 0.749 X 
10 50 erg s and T = 0.315 X 10 50 erg, respectively, at the critical rotation from 
Table 4 of lHachisu I fT986l) . Then we have n cr = 2T/J = 0.84 1 rad s" 1 , so the 
critical spin period is P sp j n cr = 2ir/Q cr = 7.47 s for a rigidly rotating 1.28 Mq 
WD. 



( > 10 9 yr) of the rapidl y rotating WD (see, e.g., Equation 10 
of lllkov & Soker1l2"0Ul) . 

4. POSSIBLE EVOLUTIONARY PATHS TO MASSIVE C+O WD 
BINARIES 

In the previous section, we present candidate binaries that 
have an evolutionary path toward a super-Chandrasekhar mass 
progenitor. Such binaries should consist of a very massive 
C+O WD (Mwd.o ~ 1.1-1.2 Mq) and an intermediate-mass 
(M2fl ~ 4-5 Mq) companion star with an orbital period of 
Porb.o ~ 0.5-1.5 days. These massive C+O WDs, however, 
have not been expected mainly because dust-driven wind 
mass-loss prevents the C +O WDs from g rowing in their mass 
in the AGB phase (e.g., iHerwig 1 120051 for a rec ent review 
of AGB stars). For example. lUmeda et al. I (Q999) calculated 
evolutions of 3-9 M Q stars for the metallicity of Z = 0.001- 
0.03 from ZAMS until the second dredge-up stage, that is, 
until the early AGB phase. They showed that an upper limit 
of C+O core mass is Mc+o = 1 -07 M Q for the initial MS mass 
of M,- = 8.1 M Q (Z = 0.02) and for M; = 7.1 M (Z = 0.004). 
If this 1.07 M upper limit is applied to all the C+O WDs, 
both the SD and DD scenarios hardly explain super-luminous 
SNe la such as SN 2007if/SN 2009dc. In this section, we 
discuss possible evolutionary paths toward our candidates. 

4.1. Metallicity Effects of C+O Core Masses 

In close binaries, the primary star fills its Roche lobe before 
it evolves to an AGB star. In relatively wide binaries, a star 
will evolve further beyond the early AGB phase (the second 
dredge-up stage). In the AGB phase, hydrogen shell-burning 
produces helium and that helium piles up on the C+O core. 
The helium intermittently ignites to trigger a shell-flash, mak- 
ing a thermal pulse phase (TP-AGB phase). The C+O core 
mass would increase further beyond 1 .07 M Q during the TP- 
AGB phase until all hydrogen envelope is lost by winds (or 
stripped away by the companion star). How much mass can 
be added to the C+O core during the TP-AGB phase ? We 
roughly assume an average increasing rate of the core mass 
from hydrogen shell-burning rate, M c ~5x 10~ 7 M Q yr" 1 , dur- 
ing the TP-AGB phase (see Equation [Q~|). On the other hand, 
the dust-driven wind mass-loss rates are observationally es- 
timated to be about M w ~ 3 x 1O" 5 M yr" 1 (Ide Beck et al. I 
2010), where the mass-loss rate levels off for stars with longer 
pulsation periods of > 850 days. Therefore, the ratio of the 
wind mass-loss rate to the mass-accretion rate is about £ = 
M w /M c ~ 60. This simply suggests that the upper limit for 
the increase in the core mass is about AM C w (My— Mi )C )/£ ~ 
6M /6O ^0.1 M Q , where M L , and M Lc are the initial mass 
and the C+O core mass of the primary, respectively. Further- 
more, recent models of massive AGB stars sho w very efficient 
third dredge-up and hot-bottom burning (e.g., Herwig 2005, 
for a recent review), both of which reduce the core growth. 
Therefore t he core of massive AGB stars seems to hardly grow 
(see, e.g., [Wsidemann 2000, for initial-final mass relation of 
single stars). 

Many authors have calculated the evolutions of TP-AGB 
stars, but full calculation through thermal pulses is extremely 
time-consuming. Thus the evolution has been often followed 
only through a small number of pulses and then replaced 
by synthetic models towards the end of the AGB. Very re- 
cently, a few full AGB calculations appeared including a 
numb er of thermal pulses until the end of the AGB phase 
(e.g.. lKarakasl2010tlsiessl2010tlVentura & D'Antonal^OOgt 
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IWeiss & Fergusori"l l2b09). The resultant WD mass, however, 
depends strongly on the adopted mass-loss law, extra mixing 
(overshooting of convection), hot bottom burning, etc. Differ- 
ent laws gave different results. 

Figure [7] shows such diverse of resulting WD masses 
through full TP-AGB calculations. All red filled symbols 
indicate the C+O core masses after the AGB evolution for 
the initial ZAMS mass of M, > 6 M . On the other hand, 
blue open symbols denote the O+Ne+Mg core masses at the 
end of the AGB evolution for the initial ZAMS mass of 
Mj > 6.5 Mq . Different s ymbols c orrespond to di fferent au- 
thors: circles dSiess 12 010). squares (Karakas 2010), tri angles 
dWeiss & Ferguson 1 120091) . and stars d Ventura & D' An tonal 
1200911201 11) . 

Ventura & D'Antona's (2009, 2011) results show ~ 
1.07 Mq as an upper limit of C+O WDs for a wide range of 
metallicities, i.e., for Z = 0.0001-0.004. On the other hand, 
Siess's (2010) results show a considerable increase in the 
C+O core mass for very low metallicities such as Z < 0.002. 
The main differences probably come from the different treat- 
ments of wind mass-loss during the TP-AGB and possibly 



extra-mixing, because TSiess I (2010) assumed Schwarzschild 
criterion of convection, which prevents the third dredge-up, 
and also assumed Vassiliadis & Wood's (1993) mass-loss 
rate, which makes the duration of the TP-AGB phase much 
longer than that of B loker's ( 1995) formalism adopted by 
IVentura& D'Antona I (12009ft . IVentura & D'Antona I ( 1201 lb 
showed that Vassiliadis & Wood's formalism gives about 9 
times more thermal pulses than Bloker's formulation for the 
same initial mass of 7 M (Z = 0.001), i.e., 348 vs. 38, be- 
cause of weaker wind mass-loss rates and thus a longer du- 
ration of the TP-AGB phase. If Vassiliadis & Wood's one is 
the case, we have C+O WDs more massive than 1 . 1 M© for 
lower metallicities, which can open an evolutionary route to 
super-Chandrasekhar mass S Ne la as massive as 2 . 3-2.7 M^. 

These t wo fo rmalisms by Vassil iadis & Wood! (119931) and 
iBlocker I d 19951) do not explicitly include metallicity effects. 
Observationally there are some suggestions on metallicity ef- 
fects of AGB winds. Local galaxies with low Z show some 
evidence for a deficienc y in the number of planetary nebu - 
lae (PNe) at [Fe/H]< -1 jMagrini et al. I2003t IZiilstra 1120041) . 
The deficiency of PNe suggests that the mass-loss rates are not 
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common envelope evolution systems with the relation of 



1 - 



>> 0.5 - 



- 



-0.5 - 




6 M e AGB -> 1.2 M. 

7 M s AGB -> 1 .2 M, 

8 M AGB -> 1 .2 M WD 



1 S WD - 
_ WD - 



2.4 



2.6 



2.8 



3.2 



FIG. 8. — Final orbital periods against the radii of AGB stars just before 
the common envelope evolution. We assume that the AGB star fills its Roche 
lobe when the C+O core grows to 1 .2 Mq and a common envelope evolution 
starts. We plot three cases of mass of the AGB star, i.e., 6 (green thin), 1 
(red thick), and 8 Mq (blue dotted). Vertical arrows indicate upper/lower 
dispersions of the factor aX in Equation (4j. (See Equations (5) and (6))- 

as high as those observed in metallicities of [Fe/H]> -1. In 
other words, a superwind do es not blow in lower metallicity 
environment (Ziilstral l2004|). These changes of b ehavior oc- 
cur at [Fe/H] f» — 1 (see also lBowen & Willsonll99ll for a the- 
oretical discussion on this metallicity effect), suggesting that, 
in lower Z ([Fe/H]< -1) environments, C+O WDs possibly 
become more massive than those in higher Z ([Fe/H]> -1) 
environments. 

From these circumstances, we expect that AGB stars pro- 
duce C+O WDs with mass higher than 1 . 1 M in lower metal- 
licities of Z < 0.002. Since the mass loss is relatively weak in 
lower metallicity AGB stars, the C+O core is likely to grow 
as massive as 1 .2 M or more before the binary undergoes a 
common envelope evolution. In such a case, the maximum 
C+O WD mass can reach 2.3-2.7 M by accretion (see Fig- 
ure^). Observational suggestions on massive C+O WDs and 
low metallicity e nviro nments of super-Mch SNe la will be dis- 
cussed in Section l531 

Finally, we comment on Nomoto & Kondo's (1991) up- 
per limit for the initial C+O WD mass: M W n.n 5. 1 -2 Mq 
for S N la explosions (see Figure 4 of iNomoto & K ondo 
[1991 . This upper limit, however, does not apply to our super- 
Chandrasekhar mass WDs. This upper limit is valid only 
when the central part of the WD stays cold so that central 
carbon burning leads to a collapse instead of explosion. In the 
normal case of Chandrasekhar mass SNe la, compressional 
heat due to mass-accretion does not reach the central part in a 
short timescale, e.g., that a WD grows from 1.3 to 1.38 Mq. 
On the other hand, for super-Chandrasekhar mass WDs, the 
timescale of mass-accretion becomes much (a few to several 
times) longer because it takes more than a few to several times 
for a WD to grow, e.g., from 1.3 to 1.8 Mq. At the time of 
central carbon ignition, the central part has been heated up, 
thus central carbon burning leads to an SN la explosion. 

4.2. First Common Envelope Evolution and Binary Orbital 

Period 

Now we examine the condition for the initial orbital period 
(forb.o = 0.5-1.5 days) to be realized when an AGB star with 
a core mass of Mq +o = 1 . 1-1 .2 Mq, unde rgoes a common en- 
velope evolution. Ide Marco et al I (1201 II) analyzed the post 



M c (M e /2+M c ) 

- G = —aO 

XR 



M c Mi (M c +M e )M 2 



2a 



f 



2a; 



(4) 



where M c and M e are the giant primary's core mass and the 
primary's envelope mass, ay and a, the final and initial sep- 
aration of the binary, M 2 the secondary's mass, which is not 
changed during the common envelope evolution, A is the num- 
ber of order unity which depends on the mass distribution of 
the primary's envelope (XR is effectively the mass-weighted 
mean radius of the envelope), and a the efficiency parameter 
of order unity. They obtained A and a for AGB stars as 



Aagb =(0.237 ±0.021) + (0.032 ±0.006) x 



M M s 

M Pl 



log a = (-1.4 ± 0.2) - (1.2 ± 0.2) log 



Mi 
Mi ' 



(5) 



(6) 



where Mms is the main-sequence mass of the AGB star, and 
Mi and M 2 are the primary and secondary masses just be- 
fore the common envelope evolution, respectively. We do 
not know the exact radius of TP- AGB stars just before the 
common envelope evolution, so we roughly assume the radius 
of R ' ~ 300-1000 Rq (see, e.g.,Figures 1 and 2 of lWebbinkl 
120071) . and the C+O core mass of 1.2 Mq, and Mi = 6, 
7, and 8 M , and M 2 = 4 M . Then we have a, ~ 700- 
2200 Rq. Substituting these values into Equations ©, 
and (O, we estimate the final separation of the binary. Fig- 
ure [H] shows the orbital period of the binary against the ra- 
dius of the AGB star just before the common envelope evo- 
lution. We obtain rather short separations and orbital periods 
of a f ~ 5-17 Rq and P orh0 ~ 0.6-3.6 days for M u = l M , 
c+o — Mwd o —1-2 Mq, and M 2 j — 4 Mq Incorporating am- 
biguity of Equations (0]i-(|6]l in common envelope evolution, 
the resultant binary periods lie possibly in the region of 0.5- 
1.5 days, which is consistent with our requirement for super- 
Chandrasekhar mass SNe la of Mwd ~ 2.3-2.7 M . 

5. DISCUSSION 

5.1. Comparison with Other SD Models 

IChen & Lil ((2009) calculated binary evolutions to SNe la 
without including mass-stripping effect. Our results can be 
compared with their cases of very sho rt WD wind phases be- 
cause of small mass-stripping effects. IChen & Lil (120091) ob- 
tained the maximum WD mass of 1 .75 Mq at t ~ 1.5 Myr (af- 
ter the companion fills its Roche lobe) for Mwd.o = 1.2 Mq, 
M 2 o = 2.5 Mq, and P or b o = 1 day (with ci = in our defini- 
tion). We obtained M WD = 1-87, 1.88, 1.89M at t ~ 1.5 Myr 
for the same initial model with ci = 1, 3, and 10, respectively. 
Our results are consistent with theirs. The small difference 
of ^0.1 Mq might stem from the difference in the mass- 
accretion rates onto the WD; Chen & Li assumed smaller 
accretion rates than ours by taking into account centrifugal 
force, thus they get a smaller WD mass. 

5.2. Detectability of SD signatures 

In a recent review, iHowell I (1201 ll) listed four major argu- 
ments strongly against the existing SD scenarios: (1) very few 
supersoft X-ray sources have been found in galaxies; (2) no 
shock interaction between ejecta and a companion star has 
been seen; (3) no hydrogen lines have been detected; (4) no 
explanation of super-Chandrasekhar mass type la supernovae 
has been given. 
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As for the first argument, although being not di- 
rectly related to the su per-Chandrasekfiar mass models, 
Gilfanov & Bo gdan I d2010t) claimed that observed soft X-ray 
fluxes of early type galaxies contradict t he SD scenari o based 
on simple assumptions. However, Hach isu et al. I (120101) 
showed that realistic evolutionary models of the SD systems 
spend a large fraction of the life time in the optically thick 
wind phase and the recurrent nova phase rather than the su- 
persoft X-ray source (SSS) phase. Thus the SSS phase lasts 
only for a few hundred thousand years. This is by a factor 
of ~ 1 shorter than those adopted by Gilfanov & Bogdan 
(2010) where the SN la progenitor WD was assumed to spend 
most of its life as a SSS. The X-ray luminosity of the SSS 
has a large uncertainty because of the uncertain atmospheric 
model of mass-accreting WDs and absorption of soft X-rays 
by the companion star's cool wind material. Thus adopting 
an average of the obse rved fluxes of existing symbiotic SSSs, 
lHachisu et al. I ((2010) showed that the observed X-ray fluxes 
in early type galaxies obtained by Gilfanov & Bo gdan I (120101) 
are rather consistent with the SD scenario. 

As for the second (and third) points, iHowell I d201 ll) 
claimed that there should be some observational evidence 
of shock in teraction between ejecta and the companion star 
( Kasenl l2010l) such as distortions in early-time light curves of 
SNe la, especially for a RG companion, but no such signa- 
tures have been observed. This is problematic for the SD sce- 
nario of WD+RG systems dHavden et al. l2010tlBianco et al. I 
20TTI see, however, a differ ent (rather contradict) result by 
Ganeshalingam et al. l (1201 lh "). These shocking signatures 
predicted in lKasenl (2010) are based on the assumption that 
the companion fills its Roche lobe. However, if the binary sep- 
aration, a, is much larger than the companion radius, R, i.e., 
a 3> R, the solid angle subtended by the companion would 
be much smalle r, a nd so would be the effe ct of shocking. 
Uustham I (120111) and lDi Stefano etal. I (120111) argued that the 
donor star in the SD scenario might shrink rapidly before the 
WD explosion, because it would exhaust its hydrogen-rich en- 
velope during a long spinning down time of the rapidly rotat- 
ing WD until the SN la explosion. In such a case, the compan- 
ion star would be much smaller than its Roche lobe, reducing 
the shocking signature, which also explains the lack of hydro- 
gen in the spectra of SNe la. 

In our super-Chandrasekhar mass SN la models, the SN la 
explosion occurs after the mass of the main-sequence (MS) 
secondary decreases to as small as 1.5-1.8 M (see, e.g., two 
evolution models in Figure [5]). Impact with such a low mass 
MS companion would not produce a large dist ortion in early - 
time SN la light curves (see, e.g., Figure 3 of K asenll2"010t) . 
so no detection of distortion is not a problem at least for our 
super-Chandrasekhar mass models. 

As for the third disfavor point, IHowell I d201 ll) claimed 
that if WDs grow to become SNe la via hydrogen accre- 
tion, there should be some indication of hydrogen in some 
form, but no such indication h as been found, e.g., in radio 
emission (Pana gia et al. 1 120061) nor i n stripped ejecta from 
the companion (Marietta et al. 1 120001) . Uustham I d2011l) and 
iDi Stefano et al. I d201 ll) argued that the donor star in the SD 
scenario would shrink rapidly to become a helium white dwarf 
before the WD explosion, having exhausted its hydrogen-rich 
envelope. Therefore, no detection of hydrogen would natu- 
rally be expected. 

In our super-Chandrasekhar mass SN la models, the SN la 
explosion occurs after the mass of the main-sequence (MS) 
secondary decreases to as small as ~ 2.0 M Q (see, e.g., model 



a in Figure [6]). Impact with such a low mass MS companion 
would lead to a strip of hydrogen-rich matter as massive as 
~ 0.1 M Q (luminous super-Chandrasekhar mass SN e la), pro- 
ducing hydrogen lines in a later phase as shown in iLeonardl 
(l2007h . 

As for the fourth (last) point, which is the theme of this 
paper, we have already shown the possible cases where the 
WD mass increases up to ~ 2.4 M Q based on our SD scenario. 
Th us, we are able to resolve the four major arguments raised 
by lHoweiil (1201 ll) against the SD scenario. 

5.3. Young Progenitors in Low Metallicity Environments 

Our calculations show that super-Chandrasekhar mass 
SNe la come preferentially from a pair of a massive C+O WD 
(> 1 Mq) and a 4-5 M main-sequence star, thus indicat- 
ing a rather y oung er population than several hundred Myr or 
so. In Section l4~lT we further show that very luminous super- 
Chandrasekhar mass SNe la such as SN 2007if/SN 2009dc 
are preferentially born in low metallicity environment because 
more massive initial C+O WDs (~ 1 . 1-1 .2 Mq) are required. 
Here we discuss these conditions in view of observational as- 
pects. 

Recent observation of the SN 2007if host galaxy suggests a 
low metallicity environment of Z ~ Z Q /9-Z Q /5 as well as a 
young population of the progenitor from the host galaxy turn- 
off mass of M ~ 4.6^j M Q dChildress et al. 11201 ll) . which 
is consistent with our luminous super-Chandrasekhar mass 
model that can attain the maximum WD mass at 4-5 M Q 
in our binary evolutions (Figure HI). We could have super- 
Chandrasekhar mass m odels even in such a lo w metallicity 
environment. Although Koba vashi et al. I d 19981) showed that 
SNe la may not frequently occur in low metallicity environ- 
ments (say, Z < Zq/10) because optically thick winds are 
weak, we here point out that the metallicity of the SN 2007if 
host galaxy is still large enough for massive WDs (with the 
initial mass of Mwd.o = 1.1-1.2 Mq) to blow strong winds 
(see Figure Q]). Even for the metallicity of Population II 
(Z = 0.001), massive WDs (M WD , > 1.1 M Q ) can blow strong 
winds and, as a result, the stripping effect works. Thus, 
in such a low metallicity environment, we expect (super- 
Chandrasekhar mass) SN la explosions in our SD scenario. 

Young populat ion of the progenitor was also suggested 
for SN 2009dc. iTaubenberger et al. I d2011l) suggested that 
UGC 10063 (suspected host galaxy of SN 2009dc if the su- 
pernova is in the edge of the tidal tail toward the galaxy 
UGC 10064) underwent a strong burst of star formation a 
few hundred Myr before the ex plosion of SN 2009dc (see also 
Figure 14 of Si lverman et al. Il201 ll for images and position 
of UGC 10063, UGC 10064, and SN 2009dc). Our favorite 
model has an MS companion of mass ~ 4-5 M . A 4-5 M Q 
MS evolves to the end of central hydrogen burning a few hun- 
dred Myr after its born. Thus our model is consistent with 
these circumstances. 

From these observational aspects, we may propose that lu- 
minous super-Chandrasekhar mass SNe la come preferen- 
tially from a low metallicity environment. This is because 
the AGB wind (or superwind) is relatively weak for Z < 
0.002 ~ Zq/ 10 so that the C+O WD can grow up to - 1 .2 M Q 
before the first common envelope evolution. On the other 
hand, the WD accretion (optically thick) wind is still strong 
enough for a ~ 1 .2 M Q WD, if the metallicity is larger than 
Z > 0.0002 ~ Z Q / 100. Considering these rare circumstances 
(a narrow range of WD masses and a low metallicity environ- 
ment), we expect rather rare events of super-Chandrasekhar 
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mass SNe la. 



6. CONCLUSIONS 



Recent observations of Type la supernovae (SNe la) sug- 
gest that some of the progenitor white dwarfs (WDs) had 
masses up to 2.4-2.8 Mq, which are highly exceeding the 
Chandrasekhar mass limit. Introducing three binary evolu- 
tion processes into our binary evolution code, i.e., optically 
thick winds from mass-accreting WDs, mass-stripping from 
the binary companion star by the WD winds, and WDs being 
supported by differential rotation, we have calculated new SD 
evolutionary models of SN la progenitors and found follow- 
ing main results: 

1. Since the mass-stripping attenuates the mass transfer 
from the companion to the WD, which prevents a binary 
even with a massive main-sequence companion of ~ 4-5 Mq 
from being merged. As a result, the WD mass can increase 
by accretion up to 2.3 (2.7) Mq from the initial value of 
1.1 (1.2) Mq. The newly obtained results satisfy the require- 
ments of recent high luminosity SNe la such as SN 2007if and 
SN 2009dc. 

2. There are three characteristic mass ranges of explod- 
ing WDs depending on the secular instability and rotation 
law: (i) Mwd > 2.4 Mq triggered by a secular instability at 
T/\W\= 0.14, (ii) M WD = 1 .5-2.4 Mq supported by differen- 
tial rotation, and (iii) Mwd = 1.38-1.5 Mq supported by rigid 
rotation. These three cases are explained as follows: 

3. If the initial C+O WD is as massive as ~ 1.1-1.2 M Q , 
r/|W| exceeds 0.14 when the WD mass reaches Mwd ~ 



2.4 Mq. We expect an SN la explosion soon after a secu- 
lar instability sets in at Mwd ~ 2.4 Mq. We regard that this 
kind of binary evolution corresponds to very luminous super- 
Chandrasekhar mass SNe la such as SN 2007if/SN 2009dc. 
These SNe la must be young (younger than a few hundred 
Myr). 

4. If the initial C+O WD mass is - 0.9-1.1 Mq, T/\W\ 
does not exceed 0.14 before Mi becomes smaller than Mb = 
3 x 10~ 7 Mq yr -1 , where we stop our binary evolution. Loss 
or redistribution of angular momentum leads to an increase 
in the central density until carbon is ignited to induce the 
WD explosion as an SN la. This waiting time for an SN la 
depends entirely on the timescale of angular momentum 
loss/redistribution. It is shorter than 10 8 -10 9 yr (angular mo- 
mentum redistribution by the Eddington-Sweet circulation in 
WDs). Therefore we suggest that these SNe la are relatively 
young, and correspond to the "prompt" component. 

5. If the initial C+O WD mass is ~ 0.7-0.8 M Q , it reaches 
Mwd ~ 1.38-1.5 Mq, which can be supported by rigid rota- 
tion. 
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